J-pairing interaction, number of states, and nine-j sum rules of four identical particles 
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In this paper we study J-pairing Hamiltonian and find that the sum of eigenvalues of spin I states 
equals sum of norm matrix elements within the pair basis for four identical particles such as four 
fermions in a single- j shell or four bosons with spin I. We relate number of states to sum rules of 
nine-j coefficients. We obtained sum rules for nine-j coefficients {(jj)J, (jj)K : I\(jj)J, (jj)K : I) 
and ((11) J, (ll)K : I\(ll)J, (ll)K : I) summing over (1) even J and even K, (2) even J and odd K, 
(3) odd J and odd K, and (4) both even and odd values for J and K, where j is a half integer and 
I is an integer. 
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I. INTRODUCTION 

The J-pairing Hamiltonian for a single- j shell is an 
important topic to study nuclear structure theory and 
also to study general many-body systems. For the case 
of J = 0, i.e., the monopole pairing interaction, the fa- 
mous seniority scheme 0,13] provides exact solutions; for 
J = •/ma, the "cluster" picture of Ref. Q presents an 
asymptotic classification of states. For other J cases, it 
was found that pairs with spin J are reasonable building 
blocks for low-lying states but little is known about exact 
eigenvalues Q . In this paper we shall go one step forward 
along this line by proving that for four identical particles 
the sum of eigenvalues for the J- pairing interaction is 
connected with a sum of nine-j coefficients. 

The enumeration of spin I states (the number of spin 
I states is denoted by .Dj in this paper) for fermions in a 
single- j shell or bosons with spin I (We use a convention 
that j is a half integer and I is an integer) is also a very 
common practice in nuclear structure theory. .Dj is usu- 
ally obtained by subtracting the number of states with 
total angular momentum projection M = I +1 from that 
with M = I Because numbers of states with differ- 
ent M's seem irregular, Dj values are usually tabulated 
in text-books, for sake of convenience. Other methods 
include Racha's method 0] in terms of seniority scheme, 
generating function mehtod proposed and studied by Ka- 
triel et al. @, Sunko and collaborators 0. All these 
works are interesting and important. However, the re- 
sults are not algebraic. It is therefore desirable to obtain 
analytical formulas of Dj. For n — 1 and 2, Di is known 
and is understood very well; but the situation becomes 
complicated when n > 3, except for a few cases with 
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Historically, the first interesting formula of Dj was 
given for the case with 1 = and n — 4 by Ginocchio 
and Haxton in Ref. . Their result was revisited by Za- 
mick and Escuderos in Ref. @. In Ref. [13, authors of 
the present paper empirically constructed Di for n = 3 
and 4, and some Z)j's for n=5. Recently, Talmi sug- 
gested a recursion formula of Dj in Ref. 11], and used 
this formula to prove D/ formulas obtained empirically 
for n — 3 in Ref. ^(|- Talmi's recursion formula is also 
readily applied to prove the empirical formulas of Ref. 
[Io| for n = 4. In Ref. we showed that Dj of n 

particle systems can be enumerated by using the reduc- 
tion from SU(n+ 1) to SO (3), and as an example, Dj for 
n — 4 was obtained analytically. 

The results of Dj for three identical particles were ap- 
plied to obtain a number of sum rules for six-j symbols in 
Appendix of Ref. Q . One can ask whether the results for 
n = 4 can be used similarly to obtain sum rules for ninc-j 
symbols. If the answer is yes, the application would be 
very interesting, because sum rules of angular momentum 
couplings is widely applied in many branches of physics, 
in particular in nuclear structure theory (angular momen- 
tum coupling-recoupling coefficients and sum rules were 
compiled in Ref. [Oj in 1988). This paper addresses the 
following question: can we obtain sum rules for ninc-j 
symbols based on studies of J-pairing Hamiltonian and 
number of spin / states for four identical particles ? Fur- 
thermore, how far can one go along this line? 

This paper is organized as follows. In Sec. II, we study 
J-pairing Hamiltonian in order to obtain summation of 
all non-zero eigenvalues in the presence of only one J- 
pairing force. In Sec. Ill, we present sum rules of nine-j 
symbols found by using these summations and number 
of states for n — 4 obtained in earlier works. In this 
paper we show that the Di formulas provide us with a 
bridge between the J-pairing interaction and sum rules 
of nine-j symbols for identical particles. The summary 
and discussion are given in Sec. IV. Appendix A present 
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formulas of nine-j symbols in some special cases. Ap- 
pendix B discusses number of matrices involved in our 
sum rule calculations. 



II. J-PARING INTERACTION 

In this section we discuss the J-pairing interaction only 
for identical fermions in a single- j shell. Similar results 
are readily obtained for bosons with spin I. Our J-pairing 
Hamiltonian Hj is defined as follows. 
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means coupled to angular momentum J and 
projection M. We take Gj = 1 in this paper. 

For n = 3, it was shown in Ref. Q that there is only 
one non-zero eigenvalue for Hj when / > j — 1, and 
all eigenvalues are zero when / < j — 1. For n = 4 
the situation is more complicated, because there can be 
many non-zero eigenvalues of spin I states for Hj, and 
most of these eigenvalues are not known except 1 = 
and / ~ /max- However, their summation is the trace 
of Hj matrix with total spin /, and is a constant with 
respect to any linear transformation. This trace can be 
obtained by summing the diagonal matrix elements 
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over K. Here J and K take only even values. This fact 
can be proved by using two-body coefficients of fractional 
parentages which are defined by 
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The trace can be calculated as follows, 
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where Gj = 1 is used. This is just the summation of 
Eq. J5J over even K. One can also regard Eq. as a 
simple generalization of the result in Ref. [3j , where it 
was shown that the non-zero eigenvalue of Hj for spin 
I states of three particles is given by the norm (j(j 2 )J : 
I\j(j 2 )J : /). Note that similar results are applicable to 
bosons with spin /. 

Let us look at nine-j symbols of identical particles un- 
der certain conditions. Based on Eq. J2J we easily find 
the following well-known fact 



(3) 



for odd I when K ^ J or K' ^ J, based on the fact that 
two identical pairs produce only even values of I and thus 
the norm in Eq. J2J) equals zero. Here K and K' take 
even values, or odd values simultaneously. This can be 
also seen from the permutation symmetry of the nine- 
j symbol, which requires the left hand side of Eq. © 
equals zero unless K + K' + 1 is even. We note without 
details that this formula is also applicable to four bosons 
with spin I, i.e., one can replace j by I in formula |J3J. 
This is a generalization of the result in Ref. 0] , where 
it was found 
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The norm of Eq. equals zero when I = 4j — 7, 4j — 
5, 4j — 4, because there are no such states. We find thus 



(4) 
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for I = 4j — 7, Aj — 5, Aj — A and J ^ K ( J, K are even) 
This is also a generalization of a formula in Ref. [l4| : 
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for I = Aj — 7, 4j — 5, 4j — 4. Similarly, we have 
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for / = Aj — 2,Aj — A. This formula was also obtained 
in Ref. and holds for both integer and half-integer 
value of j. In Appendix A, we present some explicit 
formulas of nine-j symbols with J = K = 2j or J = 
K = 2j — 1. For completeness, we also refer Refs. [4 HT3I 
IT4I ITfij , concerning formulas of six- j and nine-j symbols 
for identical particles. 

Now we enumerate the number of matrices of Eq. (0) 
with different J. This is related to the number of non- 
zero two-body coefficients of fractional parentage which 
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was obtained for specific examples in studying regular- 
ities of energy centroids in the presence of random in- 
teractions [T3. Without going to details we present the 
results of the number of matrices involved in Eq. J2J with 
different J as below. 

For I > 2j, the number of matrices with K — J is 
given by 

[(4j + 2 - J)/4] (6) 

and the number of matrices with K 7^ J is 

[(4j - J)/2] ([(4j - I)/2] + 1) /2-[(4j + 2 - 7)/4] . (7) 

The [ ] in this paper means to take the largest integer 
not exceeding the value inside. 

For I < 2j — 1, the number of matrices with J = K 
is always 1 for even 7, the case with J ^ K \s more 
complicated and the number of such matrices is given in 
the Appendix B. ft is noted that J and K take only even 
values in this Section. 

III. SUM RULES FOR NINE-J SYMBOLS 

The procedure to obtain sum rules of nine-j symbols 
in this paper is straightforward. In Sec. II we obtain 
summation of eigenvalues for H = Hj. From the sum 
rule of two-body coefficients of fractional parentage, one 



obtains 
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multiplied by the number of spin / states, 
Dj, if one sums Eq. J5J) over even J and even K, namely: 
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where Dj formulas were given in Refs. nam. New sum 

rules of nine-j symbols now can be obtained by using Dj 
formulas, Eq. J5J and Eq. 

For realistic systems both J and K are even, as in 
Eq. ((21 of Sec. II and Eq. In this paper we also dis- 
cuss sum rules of nine-j symbols under other conditions 
for J and K, such as odd J and odd K, etc. We denote 
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for sake of simplicity. The condition X of the sum rules 
for J and K will be one of the following: (1) even J and 
even K (realistic); (2) even J and odd K or odd J and 
even K; (3) odd J and odd K; and (4) both even and odd 
values for J and K. Conditions (2-4) are not physical for 
identical particles in quantum mechanics. We similarly 
define Si(l 4 , condition X on J and K) for I. 



First we present our results of ^/(j 4 , requirement X 
on J and K), which provides us with rich sum rules of 
nine-j symbols. 

For I > 2j, one obtains 

Si(j 4 ,even J even K) 
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based on Eqs. (0), ©,0, and JHJ). Let us introduce 
Io by the relation / = I max — 2Iq for even I and I = 
7 max - 2I - 3 for odd I, where J raax = 4j - 6. Using 
the Jo we can rewrite Di = D j max -2Z f° r even J and 
£>/ = -Dj max -27 -3 f° r °dd J- According to Ref. [l2|. 
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One thus has 
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The behavior of the right hand side is not easy to see due 
to those (Jo mod 6), S and [-^-j , and so on. The situation 
becomes much more transparent when one writes ^/(j 4 , 
even J even K) values explicitly. For I =even, we find 
the following formulas: 



Sj (j 4 , even J even K) = < 



for I = J max , 
for I = J nlRX - 2, 



max 

for I = J max - 4, 
for I = J max - 6, 
8 for I = J max - 8, 
10 for I = 7 max - 10, 



(12) 

for I =odd we can use Eq. I|12l) to obtain the sum rules: 
Si(j 4 , even J even K) = Si + z(j 4 , even J even if). We 
find that Si(j 4 , even J even K) has a modular behavior: 

Si(j , even J even K) 

= <%/m a x-/) mod i2)0' 4 ,even J even if) + 6 
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For I = 7 max — 1, one obtains Si (j 4 , even JK) = A based 
on the right hand side of Eqs. J5J and ©■ 

For I < 2j — 1, Eq. JSJ is less transparent to be simpli- 
fied 0], due to the complexity for Dj formula (See Eq. 
(3) of Ref. an d number of J = if and J ^ if ma- 

trices of Eq. (J2J. However, by using Eq.© of this paper, 
Eq. (3) of Ref. [T^j, and results in Appendix B, one is 
able to obtain explicitly the sum rules for I < 2j — 1; 
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has a modular behavior: 

57 (j 4 , even J even 
= % mod 12) (j 4 , even J even 7f ) + 6 
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In Eq-lHUl, m = ( j + 3/2 ) mod 3. 

In Eqs.© and ll(JH5fl J and 7f take only even values. 
It is interesting to discuss whether there are simple sum 
rules in which J and if can be both even and odd. For 
this case < 7 < 7 max = Aj. Starting from Eq. (9.29) 
of Ref. Q for J x = J 2 = J3 = J4, J12 = J13 = J, 
•hi = J2A = K, J = 7, we multiply 4(2 J + 1)(2K + 1) 
and sum over all JK (i.e., J and K take both even and 
odd non- negative integers). Using Eq. (9.28) of Ref. 0, 
we find 
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where A(JKI) means that J, K and 7 satisfy the triangle 
relation of angular momentum coupling. 

If both J and K are odd values, < 7 < 7 max = Aj. In 
this case we consider fictitious (not realistic for identical 
particles) "bosons" with a half integer spin j. According 
to Ref. [l^l j the number of states Dj for four bosons with 
spin j equals that for four fermions in a single-Z shell with 
21 = 2j + 3. As Di for four fermions in a single-/ shell 



was given in Ref. we can derive Si(j 4 , odd J odd K) 
by using Eqs. (0) and JSJl, together with Eqs. (3-5) in 
Ref. 3. Similar to Eqs. 10,(1111 and (HU), we obtain 
that when I < 2j, 
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has a modular behavior: 

SVC?' 4 , odd J odd K) 
= S (I mod 12) (j 4 , odd J odd TO + 6 

where m = (j — 3/2) mod 3 in Eq. (|17fl : and when 7 > 2j 
S7(j 4 ,odd J odd if) 
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for even i, and Si(j 4 , odd J odd if) = Si + 3(j 4 , odd J 
odd if) for odd I. For 7 = 4j - 1 (odd i), Si(j 4 , odd J 
odd if) = 0. 

For even J and odd if or for odd J and even if , < 
I < /max = 4j" — 1. For this case 

S7(j 4 , even J odd if) = Si(J 4 , odd J even if) 

= (S7(j 4 , both even and odd values for J and if) 

- Si(j 4 , even J even if) - 5/(j 4 , odd J odd if)) /2 . 

Using this relation and above results we find that when 
/ < 2j, 
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Similarly, we obtain sum rules by replacing the half 
integer j to the integer Z. First, let us study the case for 
even values of J and 7C. We find that when I < 21, 
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has a modular behavior: 
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where m = I mod 3 in Ea. H23f) . and when I > 21, 
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where 7 max = 41. 

If J and Tf take both even and odd values, similar to 
the process of obtaining Eq. (|16fl . we find for I < 21 
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We note this sum rule has the same form as Eq. I|16fl for 
I > 2j. 

For odd J and odd K values, < 7 < 7 max — 41 — 2. 
We find that when 7 < 21, 
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has a modular behavior: 

Si(l 4 , odd J odd X) 
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where m = I mod 3 in Eq. HU; and when I>21, 
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has a modular behavior 
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For odd I > 21, Sj = S I+3 with S7 max -i = 0. 

If we take odd J and even K values or we take even J 
and odd K values, < I < 7 max = 4/ — 1. For this case 

5/(/ 4 , even J odd K) = Si(l A , odd J even K) 
= (5/(/ 4 , both even and odd values for J and K) 
- 5/(/ 4 , even J even K) - 5/(/ 4 , odd J odd X)) /2 
2 [£±2] for J < 2/ , 



for I > 21 
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IV. SUMMARY AND DISCUSSION 

To summarize, in this paper we first show that the sum 
of eigenvalues of spin / states for J-pairing interaction is 
given by 

( 3 3 J 

V(l + (-) / M- 4 E( 2J+1 )( 2X + 1 )S 3 3 K 
k k { J K I 

for fermions, and 



for bosons. Then we relate them with number of spin I 
states to obtain nine-j sum rules. We study 



4(2J+l)(2if + l 



and 





4(2J+l)(2ir + l 



summing over J and K under following situations: (1) 
all J and K are even; (2) J and K can be both even and 
odd; (3) all J and K are odd; (4) J is even and K is odd. 
We also obtain formulas for special J, K and I values, 
based on the physical meaning of the norm in Eq. . 

Sum rules in Eqs. (A1-A2) of Ref. [3| can be obtained 
as a special case of the results in this paper: 1 = for 
Eqs. (|14|l and H23J) . This work is therefore a general- 
ization of some of our earlier results. We use J pairing 
interaction as a tool to obtain the sum rules but these 
results are independent of the interaction. 

In Ref. ! it was found that number of spin I states 
Di for four bosons with spin / and that for four fermions 
in a single-j shell are the same when 21 = 2j — 3. This 
produce the same value of the right hand side in Eq. (JSJ 
for fermions and bosons. Unfortunately, number of J 
and K for these two cases are different (number of J for 
bosons is l + l=j — 1/2 while that for fermions is j + l/2), 
which present different sum rules of the case with even 
values for both J and K. 

One may ask how far one can go along this line, i.e., 
to construct sum rules of angular momentum coupling 
by using formulas of £>/. As n increases, Dj formulas 
and sum of eigenvalues of spin I states become more and 
more complicated. The situation is already complicated 
for n = 4. For n = 5 there are Di formulas for only 7 ~ 
or ~ 7 max . Therefore, it is difficult to obtain Dj formulas 
and new sum rules of angular momentum couplings in 
which more particles (n > 5) are involved, except for a 
few cases with 7 ~ 7 m ay where the Dj is given by a fixed 
number series floL fll| . 
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Appendix A Formulas of special nine-j 
symbols 

In this Appendix we present formulas for nine-j symbol 

(34) 




where J = 2j or 2j — 1, based on its expansion in terms of 
six-j symbols. Value of j in this Appendix can be either 
a half integer or an integer. One sees that the nine-j 
symbol of Eq. I|34|) equals zero if / is odd, because there 
appears (— pj+4J+i — (— y phase factor if one exchanges 
the first and the second row in Ea. (|34[l . From this one 
obtains that the nine-j symbol of Eq. (|34|l vanishes unless 
/ is even. Below we discuss nine-j symbols of Eq. (|3~4~|) . 
with I being even and J = 2j or 2j — 1. 
We define 



f 

III 



3 3 2j 
3 3 2j 
2j 2j 4j — m 



(35) 



(4j + 1)2 ' 
-1 



and obtain following formulas: 
Jo _ 

n- 

fe = 
fs = 
fio 



2(4j + l) 2 (4j-l) ' 

3(2j - 1) 
2(16j 2 -l) 2 (4j-3) ' 

-3 x 5(2j - 2) 
4(4j-5)(4j-3)(4j-l) 2 (4j + l) 2 ' 

3 x 5 x 7(2j - 2)(2j - 3) 
4(4j - 7)(4j - 5)(4j - 3) 2 (4j - l) 2 (4j + l) 2 ' 
-3x5x7x9 



8(4j-3) 2 (4j-l) 2 (4j + l) 2 
(2j-4)(2j-3) 
(4j-9)(4j-7)(4j-5) ' 
3x5x7x9x11 
8(4j-5) 2 (4j-3) 2 (4j-l) 2 (4 7 + l) 2 ) 
^ (2j-3)(2j-4)(2j-5) 



We define 



(4j-ll)(4j-9)(4j-7) 



j 3 2j 
// = { 3 3 2j 
2j 2j / 



and obtain following formulas 

2 , [(2j-l)!] 2 
J0 = {-) J 7J- — TvfTT— 



(4j + l) 2 (4j-l)!2 



(2j) 



2? - [(2j-l)!f l (2j)(2j + l) 
12 1 ' (4j + l) 2 (4j-l)!2 (4j-l) ' 

= 2j [(2,7-1)!] 3 (2j)(2j + l)(2j + 2) 
74 ^ ' (4j + l) 2 (4j-l)!4 (4j-3)(4j-l) ' 



2y. 



[(2j - 1)!] 



h = (-) 



( } (4j + l) 2 (4j-l)!4 
(2j)(2j + l)(2j + 2)(2j + 3) 
(4j-5)(4j-3)(4j-l) ' 
2? . [(2j - 1)1] 7x5 



/io 



/l2 = (-) 



(4j + l) 2 (4j-l)! 16 
(2j)(2j + l)(2j + 2)(2j + 3)(2j + 4) 
(4j-7)(4i-5)(4i-3)(4i-l) ' 

2i [(2j - 1)1] 9x7 
1 ' (4i + l) 2 (4j-l)! 16 
(2j)(2j + l)---(2j + 5) 
(4i-9)(4i-7)---(4i-3)(4i-l) ' 
27 [(2J-1)!] 11 x 9 x 7/3 



(4j + l) 2 (4j - 1)! 32 
(2j)(2j + l)---(2j + 6) 
(4j-ll)(4j-9)..-(4j-3)(4j-l) ' 
2j [(2j - 1)!] 13 x 11 x 9/3 



(4j + l) 2 (4i-l)! 



:!_> 
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(2j)(2j + l)---(2j + 7) 



(4j-13)(4j-ll)-..(4i-3)(4i-l) 
fie = (-) 



.,, [(2j-l)!] 15x13x11x9/3 



(4j + l) 2 (4j - 1)! 256 
(2j)(2j + !)•■■ (2j + 8) 



(4j-15)(4j-13)---(4i-3)(4j-l) ' 

/l8 = -(-) 



, , [(2j - 1)!] 17x15x13x11/3 
(4j + l) 2 (4j - 1)! 256 
x (2j)(2j + l)---(2j + 9) 

(4j-17)(4j-15)---(4i-3)(4j-l) ' 

r/ri . 19x17x15x13x11 
J20 = (,-J . , 



5x3 



(4j + l) 2 (4j-l)! 512 
(2j)(2j + l)---(2j + 10) 
(4i-19)(4j-17)---(4i-3)(4i-l) ' 

21x19x17x15x13 



/22 = -{-Y 3 



[(2J-1)!] 



5x3 



(4j + l) 2 (4j-l)! 512 
x (2j)(2j + l)-.-(2j + ll) 

(4j-21)(4j-19)---(4j-3)(4j-l) ' 

r/o ■ 1 mi 23x21x19x17x15x13 



5x3x3 



(4j + l) 2 (4j-l)! 512 
x (2j)(2j + l)---(2j + 12) 

(4j-23)(4j-21)-..(4j-3)(4j-l) 

We define 

f J J 2j-l 

<4=< J 3 2j-l 

[ 2j - 1 2 j - 1 4j - m 

■9a = .94 = 2(4j-i)^ sec E q-( 5 ) of Sec - IL For 9m with 
larger m we obtain 



3 6 = 

g' s - 
g'w 

> 

3l2 

> 

We define 
9i 

and obtain 



3(2j-2)(16j-15) 
2(4j-5)(4j-3) 2 (4j-l) 2 
15(2j - 3)(6j - 7) 
2(4j-7)(4j-5)(4j-3) 2 (4j-l) 2 

7x5x3 
:_ (4j-5)(4j-3)(4j-l) 
(2j-4)(2j-3)(32j-45) 
4(4j-9)(4j-7)---(4j-l) ' 
9x7x5x3 



(4j-5)(4j-3)(4j-l) 
(2j + 4)(2j-5)(20j-33) 
4(4j-ll)(4j-9)---(4j-l) 



3 3 2j - 1 

3 3 2j - 1 

2j-l 2j-l I 



2 , j(4j-3)[(2j-l)!] 2 
90 [ ' (4j-l)(4j-l)! ' 



92 



( l2 , j(8i 2 -6i-3)[(2j-l)!] 2 
1 J (4i-3)(4i-l)(4j-l)! 



.94 = (-) 



2? - 3j(2j + l)(4j 2 -3j-5) [(2j-l)!]' 
(4j-5)(4j-3)(4j-l)(4j-l)! 



36 



, . 2j - j(j- + l)(2j + l)[(2j-l)!] 2 
1 ' (4j-l)! 
5(8j 2 - 6j - 21) 
X (4j-7)(4j-5)(4j-3)(4j-l) ' 
/ , 27 - j(j + l)(2j + l)(2j + 3) [(2j - l)!] 2 

■ 98 = H (4T^T)! 

35(4j 2 - 3j - 18) 

X 2(4j - 9)(4j - 7)(4j - 5)(4j - 3)(4j - 1) ' 
, , 2 , ? - j(j + l)(j + 2)(2j + l)(2j + 3) 

63(8j 2 - 6j - 55) [(2j-l)!] 2 
X 2(4j-ll)(4j-9)---(4j-l) ' 
/ ,2, 3(3 + 1)(J + 2)(2j + l)(2j + 3)(2j + 5) 

■ 912 = H W^Tji 

231(4j 2 - 3j - 39) [(2j - l)!] 2 
2(4i-13)(4j-ll).--(4j-l) ■ 

Some of above g m were also obtained for fermions in a 
single- j shell in Ref. where j is a half integer, while 
here j can be either an integer or a half integer. 



Appendix B Number of matrices with K ^ J 

for I < 2j 

Number of matrices with J = K is always 1 for an 
even value of I, which contribute 2 x ( j + ~) on the left 
hand side of Eq. (jSJ, while that (denoted by Fj here) 
with J 7^ K is rather complicated. 

For I < [j] with J > 2[(I — l)/2] and J < 2 j — 1 - 

2[f]>^ = 2[§]; 

For I < [j] with J <2[(I — l)/2], Fj = J; 
For 7 < [j] with J > 2J-1-2 [§] , Fj = [§] 4 
For [j] <I<2j with J < 2j — 1 — 2 [|] , Fj 
For [j] < I < 2j with J > 2j — 1 — 2 [7/2] and J < 

-2[7/2])l 



J-(2j-l- 



(2j-l-J)/2 + 



2[|],F J = 2i-l-2[|] 

For ft] < / < 2j with J > 2 [I] , Fj 
[1/2]. 

Because of complexity in the above classification, it is 
tedious to show ^2 JK (l + by one formula, because 

one must simplify many terms such as [ ] which means 
to take the largest integer not exceeding the value inside. 
However, one can obtain explicit sum rules by writing 
down their value and studying their individual modular 
behavior, as shown in this paper. 



